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Abstract — 

Recently, the duality between Multiple Input Multiple Output 
(MIMO) Multiple Access Channels (MAC) and MIMO Broadcast 
Channels (BC) has been established under a total power constraint. 
The same set of rates for MAC can be achieved in BC exploiting 
the MAC-BC duality formulas while preserving the total power 
constraint. In this paper, we describe the BC optimal power allo- 
cation applying this duality in a downstream x-Digital Subscriber 
Lines (xDSL) context under a total power constraint for all modems 
over all tones. Then, a new algorithm called BC-Optimal Spectrum 
Balancing (BC-OSB) is devised for a more realistic power allocation 
under per-modem total power constraints. The capacity region of 
the primal BC problem under per-modem total power constraints 
is found by the dual optimization problem for the BC under per- 
modem total power constraints which can be rewritten as a dual 
optimization problem in the MAC by means of a precoder matrix 
based on the Lagrange multipliers. We show that the duality gap 
between the two problems is zero. The multi-user power allocation 
problem has been solved for interference channels and MAC using 
the OSB algorithm. In this paper we solve the problem of multi- 
user power allocation for the BC case using the OSB algorithm 
as well and we derive a computational efficient algorithm that will 
be referred to as BC-OSB. Simulation results are provided for 
two VDSL2 scenarios: the first one with Differential-Mode (DM) 
transmission only and the second one with both DM and Phantom- 
Mode (PM) transmissions. 



I. Introduction 

In 1996, Foschini and Telatar have shown that the capacity 
of Multiple Input Multiple Output (MIMO) systems increases 
linearly with the minimum number of transmitters and re- 
ceivers [2], [3]. In multi-carrier systems with Channel State 
Information (CSI) at the transmit and receive sides, the optimal 
capacity is reached using standard waterfilling [3]. In this case, 
the preceding and equalization are given by unitary matrices 
calculated from the Singular Value Decomposition (SVD) of 
the MIMO channel for each subcarrier. In [4], this work 
has been extended to the multi-user case, where the optimal 
transmit vector covariance matrices for the MIMO MAC 
are found by a procedure called iterative waterfilling. This 
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procedure searches the optimal power allocation iteratively 
over the users using the standard waterfilling formulas over 
frequencies. 

However, finding the optimal transmit vector covariance 
matrices for the MIMO Broadcast Channels (BC) has been 
an open problem for a while. Indeed, contrary to the Single 
Input Single Output (SISO) BC where the channel is degraded 
(i.e each signal is a linear combination of the other signals 
and an additional noise), the MIMO BC is non-degraded 
which makes the capacity regions much more difficult to 
characterize. In fact, the SISO BC capacity region can be 
achieved by superposition coding at the transmit side and 
successive decoding at the receive side where the ordering 
between users is determined by the noise variance [5]. Re- 
cently, a new technique called Dirty Paper Coding (DPC) has 
been introduced [6], [7]. It can be shown that DPC achieve 
the capacity region for SISO and MIMO BC. Moreover, an 
important result has been established in the form of a duality 
theory between the MAC and BC, where the MIMO BC 
capacity regions can be characterized by their dual MIMO 
MAC capacity regions [8] . With this duality between the MAC 
and BC, it is found that the same set of rates can be obtained 
in both domains under the same total power constraint. 

However, in a multi-user scenario, and especially in xDSL, 
a constraint on the total power used by each individual 
modem is more realistic than a total power constraint for all 
modems together. In a recent paper, per-modem total power 
constraints have been applied to MIMO MAC-BC duality 
theory in a wireless context [9]. It has been shown that 
MIMO MAC-BC duality still holds if an unknown covariance 
matrix is included in the MIMO MAC optimization function. 
The capacity regions are found by means of a maximization 
on the input covariance matrices and a minimization on the 
unknown covariance matrix, requiring complex algorithms in 
order to find the optimal solution when applied to the multi- 
tone transmission. 

In this paper we devise a new algorithm called BC-OSB 
requiring less complexity compared to [9] for optimal power 
allocation under per-modem total power constraints for the 
multi-tone case (xDSL context). Our paper extends the work 
of [8] in practical scenarios where we have per-modem total 
power constraints for the BC. Indeed, it is known from [8] 
that per-modem total power constraints are not preserved by 
the MAC-BC transformations. If a single total power budget is 
used as in [8], then the power budget on each modem could be 
exceeded. Therefore we design a precoder for the BC based on 
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Fig. 1. BC (left side) and its dual MAC (right side) for tone ; 



of qij is the contribution of the user at the tone on 
the Z"' Hne. The second order moment of these data vectors 
are defined as i5[qy q^] = Qy. The covariance matrix of the 

N 

transmitted data vector _E[xiX^] = ^ Q^, since the vectors 

Qij are i.i.d. The dual MIMO Multiple Access Channel (MAC) 
for the dual uplink scenario with N users (see also Fig. 1) can 
be written as: 



Lagrange multipliers in order to transform these per-modem 
total power constraints into a virtual total power constraint and 
finally use MAC-BC transformations. This precoder rescales 
the channel matrix in order to meet per-modem total power 
constraints in the BC. The capacity region of the primal BC 
problem under per-modem total power constraints is found by 
the dual optimization problem for the BC under per-modem 
total power constraints which can be rewritten as a dual 
optimization problem in the MAC by means of a precoder 
matrix based on the Lagrange multipliers. We show that the 
duality gap between the two problems is zero. 

The paper is organized as follows. In section II, we first 
introduce the BC optimal power allocation applying this 
duality in a downstream x-Digital Subscriber Lines (xDSL) 
context under a total power constraint for all modems over 
all tones. Then, a new algorithm called BC-OSB is devised 
in section III for a more realistic power allocation under per- 
modem total power constraints. Simulation results are given 
in section IV for a VDSL2 scenario with Differential-Mode 
(DM) transmission and a VDSL2 scenario with both DM and 
Phantom-Mode (PM) transmission. 

II. BC POWER ALLOCATION UNDER A TOTAL POWER 
CONSTRAINT 

In this section, we describe a procedure for computing 
optimal transmit vector covariance matrices in the BC using 
MAC-BC duality under a total power constraint, that is under 
a single total power constraint for all tones and all modems 
in the MIMO system. This provides a reference for the 
derivations in section III. 

We assume that all modems are synchronized and use 
Discrete Multi-Tone (DMT) modulation with a cyclic prefix 
longer than the maximum delay spread of the channel. We 
consider a MIMO Broadcast Channel (BC) serving N users in 
a xDSL downstream scenario as shown on the left side of Fig. 
1. The transmission on one tone can then be modelled as: 



HfUj-l-Wi where Hf 



Yi = HiXj -I- rij Hj = 



hi 
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(1) 



where is the number of subcarriers, Xj and are respec- 
tively the transmitted data vector and the received signal vector 
of size N X 1, Hi the N X N MIMO channel matrix and 
the vector containing colored noise. The transmitted vector is 

N 

composed by N vectors of size TV x 1 such that Xj = ^ q^ . 
Qij is the vector of the j* user data where the Z* element 



(2) 

where = [un . . ■un^]'^ is the transmitted vector on tone 
i, Vi is the received signal vector of length A'', and is 
the vector containing colored noise. In this paper, we assume 
E[ninf] = I (this is without loss of generality, as correlation 
in iii cannot be exploited anyway), so that in the dual MAC 
channel E[wiwf^] = I (in the MAC, a whitening operation 
can always be applied to the received vector such that the 
noise is white). The primal problem of finding optimal transmit 
vector covariance matrices in the BC under a total power 
constraint is defined here as: 



max C^'^ 

(Qij)i=l-..Nc.J = l--.N 
N Na 

subject toY^ Trace{Qij) < P*"* 
j=ii=i 

QijhO,i = l...Nc,j ^l...N 



(3) 



where C^*^ is the weighted rate sum function for MIMO BC 

employing DPC (as DPC achieve the MIMO BC capacity 
region [8]) for a given encoding order 1,...,N-1,N (i.e. user 
1 is encoded first): 



(4) 



where ay = 1 -|-hjj( ^ Qjfc)h^. The w^ 's are the weights 

k=j+l 

assigned to the different users. In the MAC, the optimal 
detection order is actually defined by the weights and the user 
with the largest weight is decoded last [10], [11]. Assuming 
a decreasing order of weights wi > ■ • • > wk, as the MAC- 
BC duality dictates a reverse of the decoding/encoding order, 
in the BC the user with the largest weight has indeed to be 
encoded first. Thus, the first term of the sum represents the rate 
of user 1, which is encoded under the crosstalk of the other 
users. The last term of the sum represents the rate of user 
N after having removed the crosstalk from the other users. 
The weighted rate sum function is neither convex nor concave 
[8], therefore finding the optimal transmit vector covariance 
matrices in the BC is a difficult task. Fortunately, the duaUty 
between the MAC and the BC states that it is possible to 
achieve the same set of rates in both domains under the same 
total power constraint. As the optimal power allocation in the 
MAC is tractable, one can calculate optimal transmit vector 
covariance matrices in the MAC and transform these into 
optimal transmit vector covariance matrices in the BC. The 
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primal problem of finding power allocations in the MAC under 
a total power constraint P*"* is: 



(jMAC 



max 

(*i)i = l...Ne 

subject to Trace{^i) < P*' 

i=l 

*i h 0, z = 1 . . . TVc 



(5) 



with $i = E[uiuf] = diag{(l)ii, . . . the covariance 

matrix of transmitted symbols for tone i. The weighted rate 
sum function in the MAC for the decoding order N,N-1,...,1 
(i.e. user 1 is decoded last) with Successive Interference 
Cancellation (SIC) is: 



matrices found in the MAC domain into optimal transmit 
vector covariance matrices in the BC domain as follows: 



for j=N to 1 

N 

. aij = 1 + h,;,( Qifc)h,^ (OiJV = 1) 

k=j+\ 

. = I + 5: hg^ifchifc (Ba = I) 

. = B^i/^F,^Gga;/',^,^a;/'Gi^F^B: 
end for 
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MAC 



N Na 

Y.'^iY. ^°92 [det (I + hg0i,hi,B-.i)] (6) 



where B^, = I + X] hfi'/'ifchi/t- Problem (5)-(6) has been 

fc— 1 

addressed in [10], [11]. The MAC-OSB algorithm has been 
derived based on a dual decomposition approach with a 
Lagrange formulation. First, the dual objective function of (5)- 

(6) is: 



with A the Lagrange multiplier and 



max £^^^(A,(*i)i=i...;vJ (7) 



i=l 

(8) 

The dual optimization problem is: 



minunize P^^'^ (A) 

A 

subject to A > 



(9) 



By tuning the Lagrange multiplier, the total power constraint 
can be enforced. Because the dual objective function is con- 
cave with a convex constraint set, it has a unique minimum. 
As the duality gap is zero [8], [11], this minimum corresponds 
to the global optimum of the primal problem in (5)-(6). The 
search for the optimal A involves evaluations of the dual 
objective function (7), i.e. maximizations of the Lagrangian 
which, however, is decoupled over the tones for a given A. In 
particular, the Lagrangian can be rewritten as: 



with Fjj, Gij the two unitary matrices coming from the SVD 
B~^^^^h^a~^^^ = FijLy G^ and L^^ the diagonal matrix of 
singular values for user j over tone i. Note that the above 
procedure requires that the same total transmit power is used 
for the power allocation in the MAC and for the transmit vector 
covariance matrices in the BC [8]. 



Algorithm 1 Total power constraint 

1 init A = 1 

2 init step = 2 

3 init 6 = 

4 init (*i)i=i...jv„ = 

No 

5 whUe I E Trace{^i) - P*°*| > e 

i=l 

6 Exhaustive search max /^^^"^(A, (*i)j=i...jVc 

(*i)i=l--.Ne 

I if Trace{^i) - P*°* < 

i=l 

8 6=6+1 

9 A = \/step 

10 step = step - 1/2^ 

II end if 

12 A = A * step 

13 end while 

14 MAC-BC Duality 

15 for j=N to 1 Vi = 1 . . . iVc 

N 

16 aij = l + hy( E Qifc)h,^ 

17 Bi,- = I + E h^^ifchifc 

fc=l 

18 = B^V^F. .GffaV2^. .FffByV2 

19 end for 



f N r 

/:^^^(A,(*i)i=i...;vJ = E E ^jio92\ 

det (I + (/)ijhi,B->) ] - XTracei^i) J + AP 



(10) 



Therefore the maximization of the Lagrangian can be done 
by an exhaustive/iterative search on a per-tone basis [10] or 
by convex programming techniques such as interior point 
methods [11]. For further details, we refer to [10], [11]. 
Then, we can use the formulas from MAC-BC duality theory 
given in [8] to convert the optimal transmit vector covariance 



A complete algorithm description is given in Algorithm 1. 
We define e as the tolerance between the actuel total power 

constraint E Trace{^i) and the target total power constraint 

i=l 

ptot jjjg duality preserves the total power constraint, the 
transmit vector covariance matrices in the MAC (i.e. $i) or 
the transmit vector covariance matrices in the BC (i.e. Q^,) 
can be used for convergence in the loop (line 5 in Algorithm 

1 can be replaced by | E E Trace{Qij) - P*°*| > e). This 
i=ij=i 
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algorithm leads to the same rates for the different users in the 
BC (4) and in the MAC (6) domain. 



III. BC POWER ALLOCATION UNDER PER-MODEM TOTAL 
POWER CONSTRAINTS 

In the xDSL context, it is more relevant to consider a 
constraint on the transmit power of each modem separately 
instead of a constraint on the power for all modems together. 
Therefore, the goal of this section is to find optimal transmit 
vector covariance matrices in the BC under per-modem total 
power constraints, that is under a total power constraint for 
each modem over all tones. The primal problem is defined as: 



max C^^ 

(Qy )i=l...JVe,i=l...N 

N N, 



with 



subject to E E Qij.a < -P/"* V/ 
j=i i=i 

^0,z = 1...7Ve,i = l...A^ 



(11) 



with the power budget for modem I, and Qij,ii the I 
diagonal element of the transmit vector covariance matrix for 

user j over tone i coming from the /* diagonal element of the 
covariance matrix of the transmitted data vector _E[xixf^]/; = 

N 

E Qij,;(- Each diagonal element Qij,;; contributes to the 

per-modem total power constraint. C^^ is the weighted rate 
sum function as defined by (4). We aim to follow a dual 
decomposition approach similar to the approach in section II, 
and again exploit MAC-BC duality theory. The primal problem 
of finding optimal power allocations in the MAC under a per- 
modem total power constraint Pj°* is: 



max C^^C' 

(*i).=i...iv, 

N,. 

subject to E (t>,j < P*°* Vj 



i=l 

*^ >- 0,i= 1. 



2J ^ - J 

N, 



(12) 



with *j = E[viiVif] = diag{(j)i\, . . . ,4>iN) the transmit 
covariance matrix for tone i and C'^^'-' the weighted rate sum 
defined by (6). However, the MAC optimal power allocation 
computed from (12) cannot be converted directly into BC 
optimal transmit vector covariance matrices for (11) because 
the MAC-BC duality does not preserve per-modem total power 
constraints. To bypass this problem we apply a transformation 
to the dual objective function of (11) leading to an equivalent 
objective function with a total power constraint, and then we 
exploit MAC-BC duality. This transformation consists of a 
rescaling of the channel matrices by a virtual preceding matrix 
operation. The BC dual objective function corresponding to 
(11) is: 



P^'^(A)= max C^^{K{Qij)i=i...N,,j=i...N) 

(13) 



£^c^(A,(Qi,)i=i...;v„,=i 



/ N r 

-n) = E E Wjlogil 

i = 1 y J = 1 L 
Qi.h.^^"/] - Errace(AQi^)^ 



1 + 

+Trace (AdiagiPl"* , • • • > -Pl°*)) 



(14) 



with L the number of modems used in the downstream (L can 
be larger or equal than the number of active users A^)', and A 
a diagonal matrix of Lagrange multipliers diag{Xi, .. ., Al). 
The dual optimization problem is: 



minimize F^'-^{A) 
subject to A; > VZ 



(15) 



First we will solve the dual problem (15). Then we will 
show zero duality gap between the primal and the dual 
problems in (11) and (15) respectively. Rescaling the channel 
matrices by the inverse square root of the Lagrange multiplier 
matrix leads to: 




(16) 



For this equivalent channel and a given A, the dual objective 
function in the BC becomes: 



F^^(A)= max £^^(A, (Q',,.).=i...iv.,,=i...iv) 

(17) 

with 



(A, (Q^. ),=!.., 



NcJ=l 



N,, 

E 



TV r 

E Wjlog2\ 



(l + h',,Q',^.h',^a'->)] - ETraceiQ^j] 
+Trace{diag{P'T,---,P'T) 



(18) 

where a'ij = 1 + h'ij( E Q'ifc)^',^ with the rescaled 

k=j+i 

channel vectors h.' ij = h,., A'^^/^ and the re-defined trans- 
mit vector covariance matrices Q'^^ = A^/^Q^A^/^. The 
target per-modem total power constraints are also re-defined 
as diag{P'l°\ P'*/*) = A.^I'^diag{Pl°\ Pl°*)A^/\ 
One can see that (18) corresponds to (10) with A = 1, and 
so that the precoder matrix A^^/^ transforms the per-modem 
total power constraints into a virtual total power constraint by 
hiding the Lagrange multipliers into the equivalent channels 
h'ij and the new covariance matrices Q'^ . We can now invoke 
MAC-BC duality theory to transform the dual optimization 
problem in the BC into a dual optimization problem in the 
MAC for a given A under the same virtual total power 
constraint. The MAC-BC duality theory states that the sum rate 
capacity of the MIMO BC equals the sum rate capacity of the 

^For example other lines that £U"e not dedicated to specific users (or that 
are not active) are used to boost the rate of the active users through crosstalk 
coupling. 
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MIMO MAC with preservation of the total power constraint. 
Therefore F^C(A) in (13) is identical to F^^<^(A) in (22) 
resulting from the following equality concerning the sum rate 
capacity; 



matrix. We can write the equivalent primal problem in the BC 



BC 



E «V E ^^92 [ (l + h',,Q',>'^a'7.i) ] 



3 = 1 
N 



(19) 



E E ^"92 [det (l + h',^0',^.h',,.B',/) ] 



MAC 

where B'jj = I + ^'ik'P'ik^'ik- The preservation of the 
total power constraint leads to the following equality : 

N iV,. 

E E TraceiQ'.j) = E Trace{^'i) (20) 
i = ij=i i = i 

Based on the two equalities (19) and (20), the BC dual 
optimization problem in (15) can be rewritten as: 



minimize F*^'^C(A) 
subject to A( > VZ 



(21) 



with 



max £ 

(*'i)i=l...iVe 



MAC 



(A, {^'i)i=i...Nj (22) 



and where C^^'^ is the same as C^'^ in (18) except that 
we take into account (19) and (20) such that 

/ N r 

i=l \ j = l L 

det (l + h'^0',,h',,B'--i) ] - rrace(*'i))^ ^^^^ 

+Trace (dm5(PT, . . . , P'f )) 

Therefore, for a given A, we can then compute the optimal 
power allocation in the MAC by (22) and use the duality for- 
mulas of [8] to obtain the optimal transmit vector covariance 
matrices in (13). The Lagrange multipliers are then adjusted 
so that the per-modem total power constraints are enforced 
using (15). 

Up to here we found the optimal solution of the dual 
problem (15) by means of (21). In the following we will 
demonstrate that the duality gap between (11) and (15) is zero. 
Towards this end, we introduce the following theorem: 

Theorem: The primal optimization problem in the BC under 
per-modem total power constraints (11) is identical to the dual 
optimization problem in the BC under per-modem total power 
constraints (15). Hereby duality gap between (11) and (15) is 
zero. 

Proof: It is known from [9] that the capacity region of 
a BC under per-antenna power constraints for a multi-tone 
transmission is identical to the capacity region of a MAC 
under a total power constraint with an uncertain noise diagonal 



as: 

N N. 

mm 

X (*i)i=i...Ar^ J 

subject to Trace{^i) < P 

i=l 



max Y. «'•.? E ^og2 det I-hh,f<?i»jhy Yx\, 

i)i=l.. Ne 7 = li = l LV ■'/J 



(24) 

where Yx,ij = X + ^ h^^j^hifc, X is an uncertain noise 



j-i 



k=l 
:)tot 



diagonal matrix and P is the sum of the per-modem total 
power constraints in the BC. As the primal and the dual 
optimization problems of MAC under a total power constraint 
lead to the same capacity region (the duality gap is zero since 
the problem has a concave objective function under a convex 
constraint set), for a given X the dual formulation under a 
total power constraint: 



l...iVc i=l\ j = l I- ^ 



mm max 

X,A(*i 

-A Y Tracei^i) + AP*°* 

(25) 

Now, taking into account the previous definitions of the new 
covariance matrices A$i = and new total power constraint 
AP*°* = P**°*, we get: 



mm 



max g ( £ wjlog2\det(l+lhfjrijK^*£jj)] 



E rrace(**) + P 



(26) 

with Yx ij = J J2 h^0*fchifc- The interference can be 

k=l 

rewritten as: 



Y^ ., = XV2(i + 1 J2 X-V2h^0*,h,fcX-V2)xi 



/2 



(27) 



Using the following equality c!et(I -I- AB) = det(I -I- BA), 
this leads to the following optimization problem: 

Na I N r 

X,A{**)i=i...jVe i = 1 y j = 1 L 

det (l + iX- V2hg,^|.h,,X- V^B-i) ] (28) 

No \ 

- Trace(**) -|- P**°* 
i=i J 

with B;,^.. = I + i i:'x-V2h«0*,h,fcX-V2. When the 
fc=i 

uncertain noise is selected such that ^X~^/^ = A~^/^ and 

V A 

knowing that h'ij = h^A"^/^, one can easily show that 
Bx,y = B'ij with = Hence (28) will be exactly the 
same optimization problem as (21). Knowing that the (28), 
(21) and (15) optimization problems are the same, the duality 
gap between (24) and (28) is zero, and (24) has the same 
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capacity region as (11), we can conclude that the duality gap 
between (11) and (15) is zero. This completes the proof. 



A complete algorithm description is given as Algorithm 
2. We define e; as the tolerance between the actuel per- 

modem total power constraint '^[A~^^^{Y^ Q'y)A~^/^]n 

1=1 j=i 

and the target per-modem total power constraint P/"* for the 
modem. This algorithm will be referred to as BC-OSB 
(BC-Optimal Spectrum Balancing): 



Rl/R2(Mbps) 


■!U1=0.0 


wi=0.\ 


i«i=0.2 


User 1 first 


0/89.96 


136.0/89.53 


136.0/89.53 


User 2 first 


0/89.96 


136.0/89.53 


136.0/89.53 



u>i=0.3 


wi=OA 


j«i=0.5 


wi=Q.6 


136.0/89.53 


136.0/89.53 


136.0/89.53 


136.0/89.53 


136.0/89.53 


136.0/89.53 


136.0/89.53 


136.0/89.53 












u>i=0.7 


wi=O.S 


«ii=0.9 


u>i = 1.0 




136.0/89.53 


136.0/89.53 


136.0/89.53 


136.2/0 




136.0/89.53 


136.0/89.53 


136.0/89.53 


136.2/0 








TABLE II 







Rates under per-modem total power constraint in the BC. 



Algorithm 2 BC-OSB under per-modem total power con- 
straints 

1 for Z = 1 . . . L 

2 init A; = 1 

3 init stepi = 2 

4 init bi =0 

5 init (Q^j)i=i...jv,j=i...Ar = 

6 end for 

7 while 31 s.t. I E[A-i/2(E Q'v)A-i/2]n - > Q 

8 Exhaustive search max £^^^'^(A, 



9 
10 

11 

12 

13 
14 
15 

16 



(*'i)i=i.--iv, 

MAC-BC Duality 
for j=N to 1 Vz = 1 . . . TVc 



J-1 



k=j+l 



B',, =1+ Eh'^^^fch',, 
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Total power constraint 

Per-modem constraints 
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fc=i 

-1/2 



R1 (Mbps) 



100 



150 



G'„F'"B' 



17 
18 
19 
20 
21 
22 

23 end while 



end for 
f or Z = 1 . . . L 

Nc N 

if E[A-i/2(EQ%)A-i/2]n-Pr <o 

i=l j=l 
+ 1 

A; = Xi/stepi 
stepi = stept - 1/2'" 
end if 

A; = A; * stepi 
end for 
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Fig. 2. Rate region of BC-OSB in a VDSL2 system under a total power 
constraint and per-modem total power constraints. 



IV. Results 

In this section, we provide simulation results for two dif- 
ferent VDSL2 downlink scenarios. We use measured channels 
from a France Telecom binder. In the first scenario, 2 users re- 
spectively at 400 and 800 meters from the Central Office (CO) 
or the Remote Terminal (RT) are served with Differential- 
Mode (DM) lines (2x2 channel matrix). In the second sce- 
nario, 2 users both at 400 meters from the CO/RT are served 
with Differential-Mode (DM) lines and their Phantom-mode 
(PM) line (differential between the 2 common modes) giving a 
2x3 channel matrix. In this scenario, external noise is coming 
from 2 other DM lines. The spectral masks for VDSL2 Fiber- 
To-The-exchange (FTTex) are applied [12], with SNR gap 
r=0 dB (since MAC-BC duality does not hold if T >0 dB)^, 
an AWGN of -140 dBm/Hz and a maximum transmit power 
14.5 dBm per line. The frequency range is from to 



ptot 

j 



table I 

Rates under total power constraint in the BC. 



"The more general case with F >0 dB will be addressed in a future 
report. The transmit covariance matrices (Qij)i=i,,,jv^ j=i,,,Ar are optimal 
transmit covariance matrices optimal for the DPC capacity formula. Note that 
from the implementation point of view, once the optimal covariance matrices 
(Q^j)^=l...iVc,j = l...iV 'U'c determined, the transmitted data symbols can 
be constructed as follows: 
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Fig. 3. Optimal BC covariance matrix for user 1 in a VDSL2 scenario under 
per-modem total power constraints. 
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Fig. 4. Optimal BC covariance matrix for user 2 in a VDSL2 scenario under 
per-modem total power constraints 



12 MHz with 4.3125 kHz spacing between subcarriers and 4 
kHz symbol rate. The FDD band plan of VDSL2 up to 12 
MHz provides 2 frequency bands in the downlink scenario, 
namely 138kHz-3.75MHz and 5.2MHz-8.5MHz. In the second 
scenario, we transmit in the 0-30 MHz range with the FDD 
band plan of VDSL2 up to 12 MHz and transmitting in the 
all 12-30 MHz bandwidth with similar per-modem total power 
constraints. 

For TABLE I a total power constraint = 29 dBm for 
the 2 modems is chosen. The table shows the rates obtained in 
the 2x2 BC with different sets of weights (toi, w^) = (wi, 1 — 
wi). The two lines refer to two different encoding orders in 

1) The A'^ X 1 vector of the M-QAM data symbols s^j is precoded using the 
NxN Liij matrix, i.e. qjj = hijSij, such that _B[qijqg] = Qij = LjjL^^ 
from Cholesky decomposition (as Q^j is a positive semi-definite matrix). 

JV 

2) Then — ^ q^^ will be sent on the A'^ lines (the l^^ element of 

j = i 

will be sent on the l^^ line). 



the BC. The rates Rl and R2 of the two users are provided, 
with the first line corresponding to user 1 encoded first and 
the second line corresponding to user 2 encoded first. We can 
see that the larger the weight allocated to one user, the larger 
the rate allocated to this user. However, one can see that the 
difference between the possible encoding orders is negligible 
(approximately 10~^ Mbps) due to the diagonal dominance of 
the channel matrix. 

For TABLE 2 per-modem total power constraints = 
14.5 dBm are chosen. The table again shows the rates obtained 
in the 2x2 BC. As in the total power constraint case, the 
results show almost equal rates for any encoding order, due 
to the diagonal dominance of the channel matrix. Moreover, 
even varying the weights do not affect the different rates of 
the users (except for the extreme cases wi=0 and wi=l). In 
fact, the resulting transmit covariance matrices do not depend 
on the weights owing to the precoding matrix operation which 
rescales the channel vectors in the same way for these different 
weights. 

In Fig. 2, we plotted the rates of the two first tables. The 
Diagonalizing Precoder (DP) of [13] achieves 135.8 Mbps 
for the first user and 89.30 Mbps for the second user. We 
can conclude that in this scenario, owing to the diagonal 
dominance of the channel matrix, the DP achieves most of 
the capacity. Fig. 3 and Fig. 4 show the two BC optimal 
transmit vector covariance matrices as the covariance 
matrix where the tone dependency was removed. Its (m, n)* 
entry is denoted Qj.mn over all frequency tones i under per- 
modem total power constraints with m and n referring to the 
m* row and the n* column of the matrix Qij,mn- We can see 
that the optimal transmit covariance matrices have an almost 
flat power allocation in the direct channels and a much smaller 
power allocation in the crosstalk channels, where the profile 
follows the shape of the FEXT. 



Sum rate (Mbps) 


SVD TPC [14] 


SVD PMTPC [14] 


2x2 


461.46 


461.46 


2x3 


488.64 


469.97 



DP PMTPC [13] 


BC-OSB PMTPC 


448.47 


449.01 


449.01 


458.94 



TABLE III 

Comparison between different schemes 



The second set of simulation results involves a VDSL2 
scenario with a to 30 MHz bandwidth and exploiting 2 DM 
lines and their PM. In this case, the 2 DM lines are used 
for downlink transmission with external noise coming from 
2 VDSL2 DM lines of similar length, whose PSD's are set 
at -60 dBm/Hz. The PM is also used at the transmit side for 
downlink transmission, providing an overall channel matrix of 
size 2x3. By duality, this corresponds to a 3x2 MAC case 
where 2 users are transmitting in an uplink scenario and the 
2 DM lines and PM are used at the receive side. We use per- 
modem total power constraints with 14.5 dBm per modem. We 
transmit in the 0-30 MHz range with the FDD band plan of 
VDSL2 up to 12 MHz and transmitting in the all 12-30 MHz 
bandwidth with similar per-modem total power constraints. 
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TABLE 3 shows a comparison between existing algorithms 
[13], [14] for the 2 X 2 case without exploiting the PM and the 2 
X 3 case exploiting the PM. The weights are wl =w2 = 0.5. 
The algorithms SVD under a Total Power Constraint (TPC) 
and SVD under Per-Modem Total Power Constraints (PMTPC) 
provide the optimal rate sum with two-sided coordination [14]. 
Owing to the external noise, there is a rate loss between 
the SVD schemes (with two-sided coordination) and the DP 
PMTPC or the BC-OSB PMTPC (with only transmit side 
coordination). The DP and the BC-OSB algorithms lead to the 
same rates in the 2x2 case because none of these algorithms 
can properly mitigate the external noise (whitening operation 
not possible at the receive side). However, in the 2x3 case, 
the BC-OSB shows a increased rate compared to the DP. This 
is also due to the linear structure of the DP which make it 
difficult to meet the per-modem total power constraints for 
non-square matrices. 

Simulation results were performed on xDSL systems, but 
they can also apply to wireless systems. The two main 
differences between xDSL systems and wireless systems is that 
the former exhibit diagonal dominance of the channel matrix 
while the latter do not. Therefore, in xDSL systems, while 
linear precoders can achieve most of the capacity, wireless 
systems require more advanced processing algorithms (such 
as the BC-OSB algorithm presented in this paper) in order 
to achieve the maximum sum rate capacity. An extension of 
this work is to investigate the BC-OSB algorithm in wireless 
systems. 
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V. Conclusion 

In this paper we have investigated the problem of optimal 
power allocation in the MIMO BC in the context of down- 
stream xDSL. We have first described an algorithm for power 
allocation under a total power constraint, i.e a single total 
power constraint for all modems over all tones. Then, a new 
algorithm called BC-OSB algorithm has been devised for a 
more realistic power allocation under per-modem total power 
constraints, i.e a lolal power constraint I'or each modem over 
all tones, where the derivation is based on a dual problem 
formulation and an adequate transformation by a precoding 
matrix. Simulation results have been provided for different 
scenarios, namely a VDSL2 scenario with DM transmission 
and a VDSL2 scenario with DM and PM transmission. 
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